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ABSTRACT 

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m 
Shane telescope with the aim of measuring the masses of the black holes in 13 nearby (z < 0.05) Seyfert 1 
galaxies with expected masses in the range ^ 10^-10' Mq. We present here the first results from this project - 
the mass of the central black hole in Arp 151. Strong variability throughout the campaign led to an exceptionally 
clean H/? lag measurement in this object of 4.25;!;Q gg days in the observed frame. Coupled with the width of 
the H/3 emis sion Une in t he variable spectrum, we determine a black hole mass of (7.1 ± 1.2) x 10^ M©, 
assuming the lOnken et al J normalization for reverberation-based virial masses. We also find velocity-resolved 
lag information within the H/3 emission line which clearly shows infalling gas in the H/3-emitting region. 
Further detailed analysis may lead to a full model of the geometry and kinematics of broad line region gas 
around the central black hole in Arp 151. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individual (Arp 151) 



1. INTRODUCTION 



Reverberatio n mapping (iBlandford &. McKed 119821: 
lPetersonlll993h is the most successful method employed for 
measuring the central black hole mass in Type 1 active galac- 
tic nuclei (AGNs). Rather than relying on spatially-resolved 
observations, reverberation mapping resolves the influence 
of the black hole in the time domain through spectroscopic 
monitoring of changes in the continuum flux and the delayed 
response, or "echo," in the broad emission lines. The time lag 
between these changes, r, depends on the light-travel time 
across the broad-line region (BLR). Combining the radius of 
the BLR, CT, with the velocity width of the broad emission 
line gives the virial mass of the central black hole. 

To date, successful reverberation-mapping studies have 
been car ried out for app roxim ately 3 6 active galaxies (com- 
' "aLl 12004 120051) . mostly probing black 
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hole masses in the range 10' - 10^ M0. Studies of lower- 
mass AGNs have been restricted by their lower luminosities, 
requiring telescopes larger than the typical 1.5-m apertures 
that have been employed. With the goal of extending the mass 
range probed by reverberation studies, we have carried out a 
64-night spectroscopic monitoring campaign on the Lick Ob- 
servatory 3-m Shane telescope, targeting 13 AGNs with ex- 
pected black hole masses in the range ^ 10''- 10' Mq. We 
present here the first results from this project: an analysis of 
the H/3 reverberation in the nearby (z = 0.02 1 1) Seyfert galaxy 
Arp 151 (Mrk 40). Full campaign details and results will be 
presented in a series of forthcoming papers. 

2. OBSERVATIONS 
2.1. Photometry 

Broad-band Johnson B images of Arp 151 were obtained 
at the 32-inch Tenagra II telescope in Southern Arizona most 
nights between calendar dates 2008 February 26 and May 15. 
Typical exposure times were 2 x 300 s. 

The images were reduced following standard techniques. 
The flux of the AGN was measured through a circular aper- 
ture of radius 4.35", and differential photometry was obtained 
relative to 8 stars within the field. Absolute flux calibrations 
were determined on a photometric night using the Landolt 
SA-101 and SA-109 standard star fields. The calibrated fight 
curve is shown in Figure 1 . For the cross-correlation analysis, 
the B-band magnitudes were converted to fluxes. 

2.2. Spectroscopy 

Spectroscopic monitoring was carried out at the Lick Ob- 
servatory 3-m Shane telescope with the Kast dual spectro- 
graph. Arp 151 was observed from March 24 - May 20 on 
a total of 43 nights. We restricted our observations to the Kast 
red-side CCD and employed the 600 lines mm"' grating with 

A simple model of the host galaxy surface brightness profile from the 
ground-based images indicates that ~ 40% of the light within this aperture 
comes from the host galaxy starlight. 
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TABLE 1 
Light-Curve Statistics 



Time Series 


N 


■Median 


(/>" 




(1) 


(2) 


(3) 


(4) 


(5) 


Sband 


60 


L03 


1.46 ±0.04 


0.0197 


5100 A 


43 


L02 


1.58±0.18 


0.0614 


Hf3 


43 


L02 


1.14±0.19 


0.0241 



Note. — Columns are presented as follows: (1) feature; (2) number 
of observations; (3) median sampling rate in days; (4) mean flux and 
standard deviation; and (5) mean fractional error. 
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Fig. 1 . — Light curves for Arp 151. From top to bottom, the panels are: B- 
band magnitude, observed continuum flux density at rest-frame 5100 A, and 
H/3 emission-line flux. The units are as given in Table 1 for 5100 A and H/3. 

spectral coverage over the range 4300-7 100 A. Spectra were 
obtained at a fixed position angle of 90° through a 4"-wide 
slit. Exposures were usually 2 x 600 s at an airmass of ~ 1.1, 
resulting in a typical signal-to-noise ratio (S/N) of ^ 100 per 
pixel at rest-frame 5100 A in the combined spectra. The im- 
ages were reduced with IRAF and the spectra were extracted 
with a width of 13 pixels (10. 1"), with sky regions of width 6 
pixels beginning at a distance of 19 pixels. Flux calibrations 
were determined from nightly spectra of standard stars. 

To mitigate the effects of slit losses and variable seeing and 
transparency, we employed the sp ectral scaling algorithm of 
Ivan Groningen & WandersI (1 19921) to scale the total flux of 
the narrow [O III] AA4959,5007 doublet in each individual 
spectrum to match the [O III] flux in a reference spectrum 
created from the spectra obtained on the 10 nights with the 
best weather and seeing conditions. This method accounts for 
differences in the overall flux scale, as well as small wave- 
length shifts and small differences in spectral resolution, and 
has been shown to result in s pectrophotometric accuracies of 
~ 2% (Peterson et al.lll998al) . 

Finally, the spectroscopic light curves were measured. The 
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Fig. 2. — Mean and variable (rms) spectra of Arp 151 in the H/3 region. 
The naiTow lines have been removed from the rms spectrum. 

continuum flux was measured in the observed-frame window 
5187-5227 A (rest-frame 5100 A). The flux in this window is 
a combination of the AGN continuum flux, which is variable, 
and a constant component from the host-galaxy starlight. The 
H/3 flux was measured by fitting a linear continuum under 
the H/3 emission line, determined from the continuum flux 
at 4850-4900 A and 5 175-5225 A, and then integrating the 
emission-line flux above the continuum from 4900-5050 A. 
This technique includes the flux contribution from the nar- 
row H/3 emission line, which is simply a constant offset. The 
resultant light curves are presented in Figure 1, with statis- 
tical properties listed in Table 1. Figure 2 shows the mean 
and root-mean-square (rms) spectra, where the rms spectrum 
shows the standard deviation of all the spectra relative to the 
mean spectrum (i.e., the variable components of the spectra). 

3. ANALYSIS 
3.1. Time-Series Analysis 

For the time-series analysis, we place more emphasis on 
the Z?-band light curve as the driving, continuum light curve, 
although we also consider the lower S/N light curve of the 
5 100 A flux. To determine the average time lag between 
variations in the continuum flux and variations in the H/3 
emission-line flux, we follow the standard practice of cross- 
correlating the light curves using the interpolation cross- 
correlation functiori (ICCF ) method (iGaskell & Sparkel[T986l 
iGaskell & Peterson! 1 19871) as well as the discrete correlation 
function (DCF) method (Edel son & KroUklll988h . with the 
White & Peterson (1994) modifications to both. The resul- 
tant cross-correlation functions are shown in Figure 3. The 
uncertainties in the time lag are determined using the Monte 
Carlo "flux ra ndomization/random subset sa mpling" method 
described by iPeterson etal] (Il998bl |2004 . In short, the 
method samples a random subset of the data points in the 
light curves, randomizes the fluxes by applying a Gaussian 
deviation within the flux uncertainties, and cross-correlates 
the modified light curves. The procedure is carried out 1000 
times, and a distribution of lag measurements is built up. We 
include two specific measurements of the lag in Table 2: Tpeak, 
the location of the maximum of the cross-correlation func- 
tion Tmax; and Tcent, the ccntroid of the points near the peak 
of the function with r > 0.8rn,ax- The uncertainties on Tpeak 
and Tcent are set such that 15.87% of the Monte Carlo real- 
izations fall below the range indicated by the uncertainties, 
and 15.87% fall above this range (i.e., Icr uncertainties for 
a Gaussian distribution). We measure an average observed- 
frame lag of Tcent = 4.25^o gg days between the fi-band and H/3 
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Fig. 3. — Cross-coiTelation functions for Arp 151. The heavy line is the 
cross-coiTelation of the B-band flux and H/3 emission, and the data points 
show the discrete correlation between the same. The cross-correlation of 
the 5 100 A flux and H/3 (thin solid line), while noisier, is consistent. The 
dashed line is the B-band auto-coiTelation function, which peaks at zero lag 
as expected. 

light curves. The values listed in Table 2 are corrected for 
time dilation effects. 

3.2. Line Width Measurement 

The width of the broad H(3 emission line was measured in 
the mean and rms spectra. We report here two separate mea- 
sures of the line width: the full-width at half-maximum flux 
(FWHM) and the line dispersion, crime, which is the secon d 
moment of the emission-line profile ( Peterson et al.l l2004b . 
The uncertainties in the line widths are again set using Monte 
Carlo random subset sampling methods. In this case, a ran- 
dom subset of the spectra is chosen and a mean and rms spec- 
trum are created, from which the FWHM and criine are mea- 
sured. A distribution of line-width measurements is built up 
through 1000 realizations, from which we take the mean and 
the standard deviation to be the line width and its typical un- 
certainty, respectively. Additional systematic errors, such as 
those due to the exact determination of the continuum contri- 
bution, are not included in these estimates of the uncertainty. 
The line widths presented in Table 2 have bee n corrected for 
the re solution of the spectrograph following iPeterson et alj 

JIool . 

3.3. Black Hole Mass 

Following the usual assumption that the BLR kinematics 
are gravitationally driven, the black hole mass is determined 
via the virial equation 



BH 



(1) 



where r is the mean time delay for the region of interest (here, 
the H/3-emitting region), v is the velocity of gas in that region, 
c is the speed of light, G is the gravitational constant, and / 
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Fig. 4. — Velocity-resolved H/3 lag measurements for Arp 151. The error 
bars in the lag direction are the ia uncertainties on the lag measurements, 
while in the velocity direction they show the size of the velocity bin (which 
extend to the line integration limits on both sides). A clear signature of 
blueshifted gas with long lags and redshifted gas with short lags is detected, 
indicating radial infall of BLR gas. 

is a scaling factor of order unity that depends on the detailed 
geometry and kinematics of the region. 

With the premise that the Mbh - cr* relationship for lo- 
cal, c|ul^sceiU_galax^ holds for AGNs and their host galax- 
ies, [Onkenetal] (|2004) find that (/) « 5.5 for reverberation- 
based masses. This particular scaling is appropriate when Teem 
and (Jiine.rms are used for the lag and line width in the black 
hole ma ss determina tion. For the measurements presented 
here, the lOnken et al.l normalization gives Mbh = (7.1 ± 1.2) x 
10^ Mq. Individual reverberation masses, ho wever, are sub- 
ject to a typical factor of 2 - 3 uncertainty dOnken e t al.), 
likely due to differences in the intrinsic, but unknown, / 
value for each individual system. Also listed in Table 2 is 
the "virial product," assuming / = 1 . The black hole mass 
in Arp 151 is smaller than the estimate based on the stel 
lar ve locity dispersion (cr* = 124 ± 12 km s"'; Greene & H 
|2006|) and the Mbh - cr* relationship of lTremaine et alj 0200' 
which predicts 2.0^} y x 10^ M©, but the two are consistent 
within the known scatter for reverberation-based masses and 
the Mbh - o-* relationship. 

4. VELOCITY-RESOLVED TIME DELAY MEASUREMENTS 

A key goal of reverberation mapping is to recover the full 
transfer function (time delay vs. velocity structure) responsi- 
ble for the shape of the emission-line light curve in response 
to the driving continuum light curve. Determining the transfer 
function is the most promising method to potentially provide 
detailed information on the geometry and kinematics of the 
BLR. While some hints of the transfer function shape have 
been se en in certain high-quality reverberation data sets (i.e., 
iHorne et al. 1991), a fuU recovery of the transfer function has 
not yet been achieved. 

We carried out an initial analysis of the velocity-resolved 
time lag information for Arp 151 by binning the H/3 emis- 
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TABLE 2 

Rest-Frame Reverberation Measurements 



Measurement 


Value 


Units 




4 J7+0.67 


days 




■'■°'-0.24 


days 


o-lme (mean) 


1738±10 


km s~' 


FWHM (mean) 


2842 ±59 


km s~' 


ciine (rms) 


1261 ±37 


km s~' 


FWHM (rms) 


2283 ± 143 


km s~' 


CTV^/G 


1.3±0.2 


10* M0 


Mbh 


7.1 ±1.2 


10« Mq 



sion line in velocity space, where each of the eight bins con- 
tained an equal amount of variable flux. Eight light curves 
were created, one for each velocity bin, and the light curves 
were each cross-correlated with the B-band light curve using 
the methods in §3.1. Figure 4 shows the results of this anal- 
ysis: there is a clear gradient in the gas response where the 
blueshifted H/3 emission lags the response in the redshifted 
H/3 emission. This is the typical signature of radial infall: 
the gas on the far side of the AGN is moving toward us, and 
the gas on the near side is moving away from us. Outflow- 
ing (i.e., wind-driven) gas would produce the opposite effect 
(short lags blueshifted and long lags redshifted) while pure 
rotation would produce a symmetric pattern around zero ve- 
locity. Only outflow specifically precludes a determination 
of Mbh due to the non-gravitational motion of the BLR gas. 
While the kinematics of the BLR H/3-emitting gas in Arp 151 
show a strong signature of radial infall, a full two-dimensional 
echo ma p of the velo city and time delay structure in the spec- 
tra (e.g.. lWelsh & Ho rne 1991) must be recovered before we 
can fully explain the details of the BLR geometry and kine- 
matics. Such an effort is beyond the scope of this paper, but 
is cuiTently being pursued. 

The results presented here demonstrate the clearest signa- 
ture of gravitational infall in the BLR of an AGN to date. 
Some indications of infalling gas have been seen in the C IV 
broad line response in other objects, such as NGC 5 548 
(e.g.,|Crenshaw & Blackwefl 1990; Don e & Krolikll996h and 
Fairall 9 (Koratkar & GaskelL1989.) . As such, AGN BLRs 
seem to commonly exhibit signatures of infalling gas (but see 
iMaoz et al.lfT99l iKoUatschnyll2003h . The / value in Equa- 
tion 1 is directly dependent on the kinematics of the BLR, 
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thereby resulting in a different / value for an AGN BLR with 
radial infall than f or a BLR with Keplerian rotation. The 
lOnken et al.l (|2004|) value of (/) « 5.5 is determined empir- 
ically and is independent of specific BLR models. Subsumed 
into the population average, (/}, are the signatures of kine- 
matic and geometric states that are common among AGNs 
with reverberation results. While it is almost certain that the 
individual details of each BLR will result in somewhat differ - 
ent intrinsic / values for every AGN (e.g.. lCollin et alJ l2006[). 
there is currently no reason to expect that the / value for 
Arp 151 is wildly discrepant from the population average, or 
to suspect the mass derived here of having uncertainties larger 
than those typically expected for reverberation masses. 

5. SUMMARY 

We have presented the first light curves and reverberation 
analysis from our AGN monitoring campaign at Lick Ob- 
servatory. We detect a clear lag in the broad H/3 emission- 
line response to changes in the continuum flux for Arp 151, 
and we pre sent a measuremen t of the black hole mass as- 
suming the lOnken et al.l (l2004l) normalization. Initial anal- 
ysis of velocity-resolved time delays in the H/3 line shows a 
strong signature of infalling gas, but further work is needed to 
map out the detailed structure and kinematics of the BLR in 
Arp 151. 

We see strong variability in other emission lines, including 
Ha, H7, and He II, the analysis of which will be included in 
future papers. In addition, we have a Hubble Space Telescope 
Cycle 17 program to image the host galaxies of the AGNs in 
this sample, allowing us to correct their spectroscopic lumi- 
nosities for starlight and apply these new results to the low end 
of the raditi s-luminosity relationship for AGNs (iBentz et al.l 
I2OO6I l2008h . which is the primary calibration for all single- 
epoch mass estimates for broad-lined AGNs. 
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